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Abstract 

DNA has a well-defined structural transition - the 
denaturation of its double-stranded form into two 
single- strands - that strongly affects its thermal 
transport properties. We show that, according to 
a widely implemented model for DNA denatura- 
tion, one can engineer DNA "heattronic" devices 
that have a rapidly increasing thermal conductance 
over a narrow temperature range across the de- 
naturation transition (~ 350K). The origin of the 
switching behavior is the release of the base pairs 
from their confining potential as DNA denatures, 
which both softens the lattice and suppresses non- 
linear effects, increasing the conductance. Most 
importantly, we demonstrate that DNA nanoj unc- 
tions have a broad range of thermal tunability due 
to varying the sequence and length, and exploit the 
underlying nonlinear behavior. We discuss the role 
of disorder in the base sequence, as well as the re- 
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lation to genomic DNA. These results set the ba- 
sis for developing thermal devices out of materials 
with nonlinear structural dynamics, as well as un- 
derstanding the underlying mechanisms of DNA 
denaturation. 

Introduction 

Thermal transport in nanoscale materials and 
molecules has enormous potential in developing 
devices that manage heat in electronic and other 
sy stems J-^ For instance, thermal rectifiers 2 3 , ther- 
mal transistors 4 , tunable thermal links 5 , and ther- 
mal memory^ have been experimentally demon- 
strated. One can envision that many more such 
devices will become feasible as methods are de- 
veloped to engineer and control nonlinear effects 
in materials that transport heat. 

Nature has provided us with a versatile and di- 
verse nonlinear structure: DNA. Not only is the 
structural dynamics of DNA fundamentally in- 
teresting due to its relevance in biological pro- 
cesses, such as transcription 8 and replication^, but 
DNA is also being used in constructing functional 
nanoscale devices, such as a template for elec- 
tronic devices^ and molecular motors. 11 Thus, its 
ability to transport heat under different conditions 
is technologically important and may allow the 
"DNA template" to be exploited not just as a scaf- 
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fold but also as a functional device in itself. In ad- 
dition to theoretical predictions^, a recent exper- 
iment shows that incorporation of DNA into a de- 
vice can indeed give rise to nonlinear behavior in 
the thermal current.^ The experimental setup ex- 
amines a change in the thermal conductance from 
a combined duplex DNA and fluid conductor to 
a disordered single- stranded DNA layer, the latter 
being thermally insulating compared to the former. 
A complete theoretical reconstruction of the exper- 
imental results would thus need to delineate the 
role of DNA's intrinsic thermal conductance from 
the surrounding media, and examine a disordered 
layer of single- stranded DNA. 

In this work, however, we envision instead a sin- 
gle molecule of duplex DNA bridging two thermal 
reservoirs in a vapor environment. We demon- 
strate that in this setup one can tune the thermal 
transport properties of DNA by taking advantage 
of its function as the carrier of the genetic code 
via its sequence of the four bases - Adenine (A), 
Guanine (G), Cytosine (C), Thymine (T). The se- 
quence of bases determines both local structural 
properties that influence the thermal conductance 
and also where nonlinear effects give way to de- 
naturation. Together with the length of DNA, 
these characteristics make DNA's thermal trans- 
port properties highly tunable. We demonstrate 
the "engineering principles" behind tuning ther- 
mal transport and set the foundations for develop- 
ing novel thermal devices for applications in, e.g., 
nanoscale electronics. 



Theoretical Analysis 

Our starting point is the Peyrard-Bishop-Dauxois 
(PBD) model™ 5 !, which considers double- 
stranded DNA as a one-dimensional lattice of 
nonlinear oscillators. This model - the common 
model for the dynamics of DNA denaturation - 
captures the essential statistical features of DNA's 
structural transition and allows for the direct calcu- 
lation of non-equilibrium thermal transport prop- 
erties [Hia see [i] within the PBD model, the DNA 
is described by the Hamiltonian 





Figure 1: (a) Schematic of DNA between two 
thermal reservoirs. A hot reservoir (left, red) 
pumps heat into DNA, which can cause partial de- 
naturation, and a cold reservoir (right, blue) ab- 
sorbs heat from DNA. (b) Schematic of the PBD 
model for DNA dynamics. It represents DNA as a 
nonlinear lattice of fluctuating base pairs, i.e., the 
base to base distance in each base pair is confined 
in a non-linear (Morse) potential transverse to the 
backbone, neighboring bases couple only through 
the fluctuation of this coordinate. This does not in- 
clude the helical nature of the ordered state. The 
strongly fluctuating denatured DNA is predicted to 
have higher thermal conductance according to the 
common model for DNA denaturation. 

where each base pair (bp) of mass m n is repre- 
sented by stretching of its hydrogen bonds via the 
coordinate y n . The onsite and nearest-neighbor 
interaction potentials, V n and W n , depend on the 
sequence of bases. The potentials take on the 
form V n (y n ) = D n (e~ a " y " — l) 2 (known as the 
Morse potential), which describes hydrogen bond- 
ing and effective interactions due to the back- 
bone/environment, and 
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H = £ I ^ + Vn (yn) + W n (ynjn-l] 



W n (y n ,y n - l ) = ^(\+p n e-^" +y >^)(y n -y n - l ) 2 , 

(2) 

which describes the stacking interaction between 
neighboring base pairs. The scenario of interest 
is DNA in vapor similar to recent experiments, 
where it was shown that the PBD model still de- 
scribes the transition well 1819 . For the analytic re- 
sults, we consider a uniform stacking parameter, 
and then we address sequence dependent stack- 
ing using numerical results (see Ref. |20] , EU and 
the Supplemental Material for a detailed derivation 
(1) and numerical parameters). We designate the se- 
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quence of DNA by the series of bases in one of the 
strands and its complementary strand is implicitly 
implied. 

We will first focus on the thermal conductance 
ratio 



R 



where Kl(h) is the thermal conductance at low 
(high) temperatures. This quantity was introduced 
in Ref. [[2] as a way to characterize a thermally in- 
duced structural transition. Here, we will examine 
its sequence dependence. When calculated using a 
small temperature change, e.g., around the denat- 
uration temperature, it can play the role of an "on- 
off" ratio. We will see that adjusting the sequence 
and length of DNA can tune R while the sequence 
alone allows the transition temperature to be tuned 
within certain limits. 

One of the main principles behind the nonlinear 
behavior predicted by the PBD model is captured 
in the low (L) and high (H) temperature limits of 
Eq. (??): 



m n y 2 n 



7 +D;y 2 n + *-f(y n -y 



n—\, 



, (3) 



where \x = L,H. We note that the PBD model is 
an effective model of DNA near the denaturation 
transition. 22 However, these limiting forms that 
occur at much higher/lower temperatures give the 
appropriate physical description - within the PBD 
model - of DNA going from its double- stranded 
to single-stranded forms. With reservoirs attach 
on the end sites, the thermal conductance of an in- 
finite strand in these limits has the form 



K : 



J 



T H -T L 



rffi>Sr(Dj,tfJ,m,y,fi>) J (4) 



coeW 



which can be calculated analytically (see Refs. 
I12I23I24I and the Supplemental Material). Here, 
J is the heat current, the integration is over the 
frequencies that correspond to propagating modes 
(W), 7 characterizes the coupling strength to the 
reservoirs, and the transmission function X is de- 
termined by the structure of the lattice. 

|2|a) shows and R for an infinite strand with 
several different motifs, i.e., the basic unit cells of 
the DNA lattice. The high temperature conduc- 
tance of all the sequences is identical due to the 



uniform stacking interaction. The low tempera- 
ture conductance, however, varies tremendously as 
the motif is changed and is universally much lower 
than its high temperature counterpart. This behav- 
ior is driven by two distinct physical mechanisms. 
First, going from the low temperature form to the 
high temperature form results in the release of the 
onsite confining potential upon denaturation. This 
leads to softening of the phonon modes and con- 
sequently to the increase of the thermal conduc- 
tance, as described analytically in Ref. [121 Sec- 
ond, the introduction of different motifs introduces 
a mismatch between the different bases of the mo- 
tif (due to their different binding potentials). This 
results in a narrowing and splitting of the phonon 
bands, as shown in^b), and a subsequent reduc- 
tion in the low temperature thermal transport. In 
the extreme case of a semi-infinite poly(A) strand 
connected to a semi-infinite poly(G) strand the 
mismatch would be maximal. The phonon bands 
would have no overlap (see|2jb)) and the strand as 
the whole would have zero heat conduction within 
this model. Actual DNA, though, will have other 
contributions to heat conduction (e.g., the back- 
bone), which will lead to a non-zero conductance. 

A natural question to ask is what is the heat 
conductance of genomic or random sequences? 
Studying strands with periodic motifs helps un- 
derstand the behavior of random sequences. An 
infinite random sequence - and likely genomic se- 
quences - will look like a periodic strand with an 
extremely long motif. The allowed bandwidth of 
propagating modes will be narrowed by the large 
number of sites in the motif. Thus, we do not ex- 
pect heat to be conducted efficiently and the low- 
temperature K should be very small for random se- 
quences (compared to a uniform sequence or an al- 
ternating sequence). Genomic DNA, of course, is 
always finite and not completely random. We ex- 
pect, however, that the result will be similar to that 
of a random sequence. However, small regions of 
the genome can look very different from a random 
sequence, and nature may exploit sequence varia- 
tion to optimize heat (signal) transport. 
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Figure 2: Engineering the thermal conductance of DNA. (a) The thermal conductance ratio and con- 
ductance for infinite strands with different periodic motifs at the high/low temperature limits of the PBD 
model. The low temperature conductance (blue circles), Kj,, drops precipitously as the motif of the se- 
quence is enlarged while the high temperature conductance (red squares), Kjj, stays the same. The thermal 
conductance ratio, R, introduced in Ref. 12, is the ratio of the high to low temperature thermal conduc- 
tance and characterizes how DNA's thermal conductance changes as it denatures, (b) Calculated phonon 
bands, from the lattice parameters used to fit the PBD model to experimental denaturation curves, for 
selected DNA sequences at low temperature. The left most bands are for a homogeneous sequence of G 
(upper, purple band) and a homogeneous sequence of A (lower, green band). As the motif is enlarged, the 
bands will become more narrow and also split. In other words, periodic sequences with increasing lengths 
of periodicity will decrease the bandwidth of the phonon modes, drastically reducing DNA's ability to 
conduct heat at low temperature, (c, d) Engineering the operation temperature of a DNA-based thermal 
device using the sequence dependence of the denaturation transition, (c) Heat current across select DNA 
sequences with AT '=10 K and with 90 bp length, where the first and the last 20 bp segments are connected 
to Langevin reservoirs, (d) Heat capacity normalized per number of base pairs for a 90 base-pair (bp) 
strand with periodic boundary conditions. The inset shows the heat capacity for an A200G200 repeat unit. 
The peak in the heat capacity marks the spot where DNA transitions into the denatured regime, and thus 
it sets the temperature where the increase in conductance ends, (e, f) Engineering the range of operation 
using the length of DNA. (e) Heat conductance of poly(A) of length 50, 100, and 200 bps (with an ad- 
ditional 20 bp segment connected to a reservoir at each end), (f) Heat capacity around the denaturation 
temperature for poly (A) of varying lengths. 



Numerical Results and Discus- 
sions 

The PBD model was developed to describe the 
properties of DNA around the denaturation tran- 
sition. The analytic results above addressed high 
and low temperature limits (within the simplifi- 
cation of uniform stacking interactions). In or- 
der to understand the extent to which the values 
of R realized in [2fa) can be realized in a nar- 
row temperature range around a transition, we per- 
form numerical simulations of the full model in- 
cluding a sequence-dependent stacking interaction 
(parametrized in Ref. |2T|). The heat current is ob- 
tained by keeping the temperature difference be- 
tween the heat baths constant (7# — Ti = 10 K) 
and scanning the average temperature (T) = (7# + 
T L )/2. We consider 90 base-pair (bp) long strands 
with 20 bps at each end connected to Langevin 
reservoirs. The damping of the individual sites by 
the reservoirs is 0.5 ps _1 . This is large enough 
to keep the very ends at the temperature of the 
reservoirs, while still allowing the sites to fluctu- 
ate at their natural frequency. The Supplemental 
Material has further details on the numerical sim- 
ulations. 

|2fc) shows the heat current for several se- 
quences. The analytical results above predicted 
that the heat conductance at temperatures higher 
than the melting temperature is insensitive to a 
particular sequence since the stacking interaction 
was assumed uniform. This is not the case in the 
numerical simulations, where the stacking poten- 
tial assumes a more realistic sequence-dependent 
form. Accordingly, the heat current exhibits a 
dependence (although rather weak) on sequence 
in the temperature range 400 — 450 K. However, 
the heat conductance of DNA increases drasti- 
cally when the temperature crosses the denatura- 
tion point. That this is indeed the denaturation 
transition where the conduction of DNA strand 
changes rapidly (versus temperature) can be seen 
by the correspondence between these curves and 
the peaks in the heat capacity shown in |2jd). 
Around the transition, poly(AG) and poly(A2G2) 
have about the same conductance and heat capac- 
ity, implying that the denaturation bubbles - where 
the two strands locally come apart - are much 



longer than the motif and, thus, only the average 
sequence matters. Thus, the sequence will change 
the denaturation temperature and tune the "oper- 
ating temperature" of the device. The GC base 
pair has a higher dissociation energy than AT, and 
thus its incorporation into a strand increases the 
transition temperature. 

Furthermore, while the sequence can change 
R measured by the high/low temperature limits 
and the "operating" temperature, the on-off ratio 
around the "operating" temperature does not dras- 
tically change. However, tuning the length of the 
DNA nano-j unction allows one to tune this im- 
portant device characteristic. [2je) shows the heat 
conductance of poly(A) strands of various lengths. 
Below the transition, the strand is anharmonic and 
is expected to demonstrate finite well-defined con- 
ductivity, i.e., the conductance is expected to be 
inversely proportional to the length of the strand 
so long as this length is longer than the typical 
bubble size. As seen in[2je), near (and above) the 
transition, the conductance weakly depends on the 
length, signifying that the harmonic, high tempera- 
ture Hamiltonian is being approached. This is fur- 
ther supported by [2]T), where above the transition 
point the heat capacity is seen to rapidly approach 
C/ks = 1 - the harmonic limit. 

At temperatures further below the transition, 
the conductance drops inversely with length of 
the DNA within the error of the simulations^. 
This observation is in agreement with the heat ca- 
pacity which shows the transition narrowing for 
longer strands of DNA. This is simply an indica- 
tion that as the bubble length becomes comparable 
to the strand length, denaturation has effectively 
occurred and this smears out the transition temper- 
ature range. We conjecture that using sequence ef- 
fects (e.g., the suppression of the low temperature 
conductance shown in [2]^a)) together with length 
will allow for even more drastic on-off jumps in 
the conductance. However, a more detailed study 
of DNA, including backbone effects, will be re- 
quired to investigate this issue. 

Conclusion 

To summarize, we have examined the thermal 
transport properties of DNA as described by the 
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PBD model and predicted that DNA nano-devices 
act as a thermal switch - where the thermal con- 
ductance rises many orders of magnitude as the 
temperature of the DNA strand is driven across 
the denaturation transition. The operating prin- 
ciple behind this behavior is the release of the 
base pairs from their confining potential, which 
both softens the lattice and suppresses nonlinear 
effects.^ Using analytic calculations and numeri- 
cal simulations with sequence-specific parameters, 
we have shown that the operating temperature of 
the thermal switch can be tuned by choosing dif- 
ferent DNA motifs and that the "on/off" ratio can 
be tuned by the DNA length. 

Further possibilities for engineering thermal 
transport may be offered by molecular or chemical 
modification of the nucleotides, using much longer 
sequences (e.g., see the inset of [2jd), showing that 
a two step jump in conductance may be possible), 
and exploiting extrinsic changes in heat conduc- 
tion (e.g., due to a structural change modifying the 
surrounding environment in addition to changing 
intrinsic properties, as in a recent experiment*^). 
This work sets the foundation to developing ther- 
mal switches out of materials and molecules with 
nonlinear structural dynamics. In addition, it will 
allow one to test underlying mechanisms for struc- 
tural transitions^ and, in particular, the dyna mical 
behavior captured within the PBD model.^GZEE] 
We speculate that biological systems may take ad- 
vantage of such nonlinear behavior in engineering 
their own control of heat flow and signaling. 
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BASIC MODEL AND PARAMETERS 



II. ANALYTIC DERIVATION OF THE 
THERMAL CONDUCTANCE 



In the main text we discuss the thermal transport 
properties of a DNA nano-junction. We designate the 
sequence of DNA by the series of bases in one of the 
strands and its complementary strand is implicitly im- 
plied. Furthermore, a DNA sequence with a periodic 
motif is denoted by, e.g., poly(AG), which is the alter- 
nating sequence AG AG AG. . . 

We describe the DNA in the framework of the PBD 
model [1-3] widely used to describe the dynamics of DNA 
denaturation. Essentially, the DNA is represented as a 
set of balls and springs with the balls representing the 
base-pair stretching and the springs encoding the com- 
plex interactions of bases within pairs (H-bonding) and 
between pairs (stacking). The PBD Hamiltonian is 

,-,2 



+ V n (y n ) + W n (y„,y n -i) 



(1) 



y n is the stretching of the base pair hydrogen bonds, and 
m n are the masses of the base pairs, which are assumed 
uniform throughout the strand with m n = m = 300 
a.m.u. [4, 5]. The functions V (y n ) and W (y n , y n -i) 
describe the effective (intra-pair) potential felt by the 
hydrogen bond and the effective inter-pair stacking in- 
teraction, respectively. These potentials take the form 

V n (y n ) = D"(e-°»»» - l) 2 

K n 



W n (y n ,y n - 



-(l + Pn e-^+yn- 1 ) ){yn _ yn _ 1 ^ 



(2) 

Recently developed parameters [4, 5] of the Morse 
potential for complementary base pairs are D A ~ T = 
0.05 eV, cla-t — 4.2 A for the A — T base pair, and 
D G - C = 0.075 eV, a G _ c = 6.9 A" 1 for the G - C base 
pair. The parameters of the stacking interaction be- 
tween successive base pairs depends on both the pairs 
and their orientation. However, we first consider the 
uniform "average" stacking potential with parameters 
K n = K = 0.025 eV/A~ 2 , p n = p = 2, and (3 = 0.35 A" 1 
[5] . The effect of the sequence-dependent stacking inter- 
action is addressed in the numerical simulations using the 
Langevin equation. 



The thermal conductance of a classical harmonic lat- 
tice can be found analytically. Our starting point is 
to consider the limiting cases of the single-coordinate 
Hamiltonian (Eq. (1) in the main text) for a lattice of 
length N. The details of the high and low-temperature 
expansions were presented in Ref. [6]. The low- (L) and 
high- (H) temperature limits can be approximated by a 
harmonic Hamiltonian of the form 



m Vn , n n 2 . i v 



(3) 



where \i = L, H, and the corresponding coefficients 
K H = K, D H = and K L = K(l + p), D L = Da 2 . 

Once a harmonic Hamiltonian is obtained, one can fol- 
low the procedure of Refs. [7, 8]. To simplify the analysis, 
the lattice is coupled to two heat reservoirs at the first 
and last sites, which gives the equations of motion 

my n = -2(D™ + K^)y n + K^(y n ^ + y n+1 ) 



+ {S n ,l + S n ,N t ) 



f dt'A{t-t')y n {t')+r] n {t) 

J — OO 



(4) 



We choose the spectrum of the dissipation to be ohmic, 
A(uj) = —ijuj, with coupling 7, and the noise to be a 
white noise, (Vl/h(^)Vl/h{^')) = 4irT L/H -fS(uj + u>'), 
with T L / H the low and high reservoir temperatures. We 
consider a periodic lattice of length N t = N x N s , where 
N s is the length of the motif. This form for the reservoirs 
will satisfy the fluctuation-dissipation theorem. The re- 
sulting equations of motion are 

my n = -2(£>™ + K^)y n + ^(y n _i + y n+1 ) 

+ (<W + 5n,N t ) [-lVn{t) + Vn(t)} ■ (5) 

The solution for the coordinates has the form 

/•OO 

y n {t) = (1/27T) / dwY-^(u;)fj m (u;)e iut , (6) 



where f\ is a vector of length Nt with the first and last 
components being t]l/h{u) and the rest being zero. It 
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represents the coupling of the reservoirs to the ends of the 
lattice. The N t x N t matrix Y = (f> — lo 2 M — A encodes 
the solution. Here cf> nm = 2(D^ + K^)5 n ^ m - K fl S rhm+1 - 
^A,m-i ; Mij = mSij, and A u = A NtNt = A(w) and 
Anm — otherwise. 

The heat current flowing into the lattice is J = 
([/_ dt'A(t — t')yi(t')]yi(t)), where the average is over 
the noise. Setting 7 = Am, the heat current becomes 

J„ = m f dLJU 2 {{V liNt - A^m 2 ^,^) 2 

T J-oo 

+ \ 2 Lj 2 m 2 (V hNt _ 1 +V 2 , Nt ) 2 }- 1 \G\. N f, (7) 

where AT = Tjy — Tr is the temperature difference of 
the reservoirs, C\ t N t is the cofactor of Y^jv, and V n ^ m is 
the determinant of the submatrix of (<j> — u> 2 M) from the 
n-th row (column) to the m-th row (column). It follows 
that \C hNt \ 2 = K 2N t- 2 and V n;m = K n ~ m ^V^ m . 



The elements 



2N t 
u 2.2N t 



3^ Nt ~ X ) = %\ N , where 

U 2,2N t - 



T s = T\Ti ■ ■ ■ Tn s is the composite transfer matrix of one 
segment with length N s . Since each transfer matrix % is 
unimodular, for an infinite lattice the allowed propagat- 
ing modes correspond to T s with eigenvalues e ±lq . We 
notice that for those propagating modes, 



U l,N s 



U 2,N S 



2cos(g). 



(8) 



This equation determines uj(q). In general, for a ba- 
sis with N s elements, there are N s bands that satisfy 
Eq. (8) though the bandwidths can be substantially re- 
duced as N s increases. Moreover, [7^]^ = [cos(./V<7)]l + 
[sin(Nq)/ sm(q)][T s - (cos(g))l], where 1 is the 2 x 2 
identity matrix. Then we rewrite the denominator of 
Eq. (7) as \z A \ 2 , where z A = (1 + \ 2 u 2 m 2 / K 2 ) cos(Nq) + 
(sin(Nq) / sm(q))(z c - (1 + \ 2 u?m 2 / K 2 ) cos(q)) and z c = 
(P? iJVs - \ 2 u? m 2 V\ Ns _jKl) - iXuim/K^Vl^ + 
T>2 N )■ Eq. (8) determines uj{q) so the integration of 
Eq. (7) can be considered as an integration over q for 
those propagating modes. As N — > 00, one uses the 
formula J Q n dqF(q,Nq) ->• ± f " dq f Q n dxF(q,x) by 
treating x = Nq as an independent variable in the 
N — > 00 limit. After integrating over x and q, one obtains 
the thermal current. 

For a uniform lattice, N s 



1 and D" 



Dp so 



This gives for the low and high temperature thermal con- 



ductance, k m = J^/AT, 



k B mK 2 
47 s 



2 7 2 , 2 7 2 D p 



mK.j 



mK 2 

A* 



with 




m 2 Ki 



(11) 



(12) 



With these expressions one can explicitly find the thermal 
conductance ratio R. We have verified that for reservoirs 
contacted to a single site on each end, the thermal con- 
ductance from our numerical simulations agree with our 
analytic formula to within 10 — 15%. The error may be 
attributed to finite size effects in the numerical simula- 
tions. 

We can take various limiting forms of these equations. 
If we define the prefactor as and a dimensionless reser- 
voir coupling as 



7m 
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the expressions for the conductance become 



2 ^ B„ 



l + 2 7 2 + 2 7 2 -^ 



(13) 



(14) 



with 



D„ 



B » = V 1 + K + 4 7^ + + 4 7 



The appropriate limiting forms for our case are the fol- 
lowing. When the high temperature harmonic limit has 
no onsite potential, then the heat conductance becomes 



k h = K H 



1 + 2 7 | - y/l + bfe 



(16) 



For the low temperature limit that has a much greater 
onsite term than the nearest neighbor coupling, i.e., 
Kl/Dl <C 1, the heat conductance becomes 



D L ' 



(17) 



7; = 71 



2(1 + - (m/K^u 2 -1 



1 







(9) 



Eq. (8) gives 2cos(<z) = 2(1 + D^/K^) - (m/K^oj 2 . 
After changing variables from u> to q that satisfy this 
constraint, the final expression (for an infinite lattice 
(N -> 00)) is 



AT 



27TTO 



/•27T 

Jo ' 



dq- 



sin' 



: (<z) 



2 7 2 
mK u 



cos(q) 



which also assumes that the dimensionless coupling to the 
reservoirs is 7^ > 1. For strong coupling to the reservoirs, 
the ratio becomes 



R. 



2K H D L 



(18) 



This is the analytic expression that demonstrates that 
softening of a harmonic lattice increases the thermal con- 



(10) ductance ratio (Dl 



[). The strong coupling limit 



gives the extreme value of R. 
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Now we show that the characteristic frequency in the 
PBD model is lowered as T crosses T c from below. 
For the low temperature Hamiltonian, the correspond- 
ing equation of motion is 

my n = -{2Da 2 y n + K(l + p)[(y n -y n _ 1 ) + (y n -y n+1 )]}. 

(19) 

From the ansatz y n = y°e iwt lkn , one obtains the 
phonon spectrum as 

muj 2 = 2Da 2 + 2K(l + p) [1 - cos(fc)] . (20) 

Thus, the frequency band of phonons is y/2Da 2 /m < 
lo < yJ[2Da 2 + 4K(l + p)]/m. For the high temperature 
Hamiltonian, the equation of motion is 

my n = -K[(y n - y„_i) + (y n - y n+ \)] (21) 

and the phonon spectrum is 

muj 2 = 2K[l - cos(fc)]. (22) 

The frequency band of phonons is < u < -J 'AK/m. 
The two limiting Hamiltonians are both harmonic and 
the characteristic frequency is indeed lowered, and sim- 
ilar considerations apply for other models. In the low 
temperature limit, the onsite potential stiffens the DNA 
compared to the high temperature limit, which results in 
the raising of the phonon spectrum of the low temper- 
ature limit compared to the high temperature. In the 
latter, as well, the nearest neighbor coupling drops from 
K (1 + p) to K shrinking the bandwidth. This trade-off is 
responsible for the change in thermal conductance across 
the transition. If the drop in nearest neighbor coupling 
is small, then the softening will dominate, and the heat 
conductance will increase because these softened modes 
can conduct heat more effectively. 

Next we consider a lattice with periodic segments. For 
a periodic lattice with alternating (AG) pairs, N s = 
2 so there are two onsite parameters D 1 ^ 2 and % = 
(V% -X>° 2 \ ([ Xl ] -l\([x 2 ] -l\ 

where \x x \ = 2(1 + j^) - f-u; 2 and [x 2 ] = 2(1 + - 
j^-Lo 2 - From Eq. (8) and | cos(g) | < 1 one can find two 
bands which satisfy the condition: 

Ubi <u 2 < Ub2 and u t \ < to 2 < u t2 . (23) 

Here u bl = (K^/m)[2 + (£>* + Ify/K^ - 
^+(D^ 1 -Df l ) 2 /K 2 )}, u b2 = 2{K ll /m){\ + D^/K^, 
uti = 2{K ll /m){\ + Dl /Kp), and u t2 = 
(K,/m)[2 + (Dl + DD/K, + ^4 + (D* D^/K 2 )}. 



One can show that the current is 

4r = X % I ^|^sin( g )|[(l + AVm 2 /^) x 

\m f f Ub2 f u t2\ 

X 2 m 2 u/K 2 )(4 + IDHKp + 2£ 2 /i^ _ 
2mu/K„)}-\ (24) 

Here we have used u — uj 2 and cos(q(u)) — (1/2) [(2(1 + 
Dl/Kp) - mu/K„)\\2{\ + Tfi^K^ - mu/K^} - 1 from 
Eq. (8). 

The heat currents for more complicated bases can be 
derived in a similar way. For example, the heat current 
for a poly(A 2 G) strand is given by 

ir = ^l «Wi-«rf(ff)[(i + 

AT 2ttK m J weW 

\ 2 m 2 u/ K 2 )\[xi}[x 2 ] + [ Xl ] 2 - 2I]- 1 (25) 

Here the three conduction bands are determined by 
2cos(g) = D° 3 - X>° 2 = [xi] 2 [x 2 ] - 2[ Xl ] - [x 2 ]. The 
heat current for a poly(A 2 G2) strand is 

— = — — / dnVl-cos 2 (<z)[(l + \ 2 m 2 u/Kl) x 

\[xi] 2 [x 2 ] + [x^] 2 - 2\ Xl ] - 2[x 2 }\}- 1 (26) 

The four conduction bands are determined by 2 cos(q) = 
^14-^2 3, where 2?° 4 = [ Xl \ 2 [x 2 ] 2 - [ Xl ][x 2 ] - [xi] 2 - 
[x 2 ] 2 + l, P° 3 = [xi] [x 2 ] - 1, 3 = [xi] [x 2 ] 2 - [xi] - [x 2 ] , 
and V\ A = [xi] 2 [x 2 ] - [ Xl ] - [x 2 ]. 

For parameters relevant to real DNA, we have evalu- 
ated the thermal conductance (in the harmonic limit) 
of poly (A), poly(G) homogeneous lattices as well as 
poly(AG), P oly(A 2 G), poly(AG 2 ) and poly(A 2 G 2 ). They 
qualitatively agree with our numerical simulations. 

III. SIMULATION DETAILS 
A. Transfer matrix 

The transfer matrix formalism is based on the possi- 
bility to evaluate the partition function for a classical 
non-linear lattice using the matrix algebra [9]. Specifi- 
cally, the total classical partition function for the PBD 
model is given by 

Z = Z T Z V , (27) 

where Z T = (2irmkBT) N / 2 is the "kinetic" partition func- 
tion and the "potential" partition function reads as 

Zu = [ ft ^ e-^v«)+^v»>v»-i)W>B T . (28) 

J 71=1 
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The continuous integration can be represented through 
summation on the grid yielding for periodic boundary 
conditions (i.e., t/o = Vn) 



Ztt = Tr 



N 



II Mn 



(29) 



where matrix M n is defined through its matrix elements 
as 

[M n ]ij = Axel^ Vn{xi)+Wn{xi ' x ^ + i v "- l{x ^ /kBT , (30) 

where Xi is the coordinate on the grid with the grid step 
Ax, i.e., Xi = xq +i Ax. The parameters of the grid used 
in our calculations are x m i n = —1-5 A, x max = 50 A with 
Ax = 0.05 A. 

In the case of a homogeneous lattice, i.e., if V n (y n ) = 
V(y n ) and W„(y n ,j/„_i) ee W{y n ,y n _ 1 ), all the ma- 
trices M n are identical and symmetric, which results 
in the possibility to evaluate the trace in Eq. (29) as 
Zjj = J2 i \™ , where A* are the eigenvalues of symmetric 
matrix M n . Naturally, if the lattice becomes very long, 
only the largest eigenvalue, \ max of the transfer matrix 
contributes to the partition function yielding Zjj = X^ ax 

The natural generalization of this procedure to the case 
of DNA strand made of repeating motif is as follows. 
First, matrices corresponding to a single repeating unit 
(r) are multiplied yielding M r = Y\ ner M n . The poten- 
tial partition function is given by Tr[M^ r ], where N r is 
the number of repeat units in the strand with periodic 
boundary conditions. Therefore, the potential partition 
function becomes Zjj — . vf r , where Vi are the eigen- 
values of M r . We note here, that even though matrix 
M r might not be symmetric, the eigenvalue problem for 
such a matrix is still well defined and eigenvalues is all 
what is needed to evaluate the "potential" partition func- 
tion above. Specifically, it can be shown that the Schur 
decomposition, which yields eigenvalues, is always pos- 
sible for any square matrix. This decomposition does 
not guarantee the completeness of the eigenvector-based 
basis, but only knowledge of eigenvalues is required for 
the evaluation of the potential partition function for the 
lattice with periodic boundary conditions. 

Once the partition functions are found, the heat capac- 
ity of the DNA strand, normalized per number of sites, 
can be evaluated as 

T d 2 F 

c = -nW < 31) 

— kBT]n(Zr^u) ' 1S the free energy of the 



where F = 
DNA strand. 



B. Langevin dynamics 

To study the dynamics of the DNA out of equilibrium 
we solve numerically the Langevin equation, which de- 





K(eV/k 2 ) 


AA 


0.023 


GG 


0.019 


AG 


0.0232 


GA 


0.0185 



Table I. Harmonic parameter of the sequence-dependent 
stacking interaction of the PBD model. 



scribes the dynamics of a Hamiltonian system in the pres- 
ence of thermal baths. The Langevin equation is given 
by 



my n = - 



dW 
dy n 



dV 

dy n 



T n y n + f(t), 



(32) 



where W(y n ) and V(y n ) are the potentials described in 
Eq. (1) of the main text. The DNA strand is split into 
three regions, the two ends, each of length I, serve as the 
Langevin thermal reservoirs at temperatures Tl and Th ■ 
This means that the friction term T n only operates for 
n within the thermal reservoirs. The fluctuating term 
f(t) is Gaussian white noise which obeys the fluctuation- 
dissipation relation (f(t)f(t')) = 2Y n k B T L{H) 5{t-t') for 
the low and high temperature reservoirs, respectively. In 
our simulations, the parameters of each Langevin ther- 
mostat have always been set to I — 20 and T n /m = 
0.5 ps- 1 . 

The middle region of the length M is the free DNA 
strand, which is driven out of equilibrium by the 
Langevin reservoirs when Tl ^ Th- The length of 
the middle region has been varied in the range M = 
50 — 200 in our simulations. The parameters of the 
PBD model are given in the main text except for the 
sequence-dependent stacking interaction strength. These 
are adopted from Ref. 5 and compiled in Table I. The 
equations of motion are integrated with the fourth-order 
Runge-Kutta method. The local heat current is given by 



Jn 



9W(y n ,y n -!) 



. The simulations are performed 



long enough to allow the system to reach its steady state, 
where J n does not depend on n, as long as n is within 
the "free" middle region of DNA. 
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